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Abstract. We consider concentric circular defects in the two-dimensional Ising model, 
which are distributed according to a generalized Fredholm sequence, i. e. at exponentially 
increasing radii. This type of aperiodicity does not change the bulk critical behaviour 
but introduces a marginal extended perturbation. The critical exponent of the local 
magnetization is obtained through finite-size scaling, using a corner transfer matrix 
approach in the extreme anisotropic limit. It varies continuously with the amplitude of 
the modulation and is closely related to the magnetic exponent of the radial Hilhorst-van 
Leeuwen model. Through a conformal mapping of the system onto a strip, the gap- 
exponent relation is shown to remain valid for such an aperiodic defect. 



1. Introduction 

The influence of a layered aperiodic modulation of the couplings on the critical behaviour 
of the two-dimensional Ising model has been much studied recently (see [?] and references 
therein). Such a perturbation may be relevant, marginal or irrelevant, depending on the 
sign of the crossover exponent [?, ?] 

$ = 1 + 1/fw - 1) , (1.1) 



which involves the bulk correlation length exponent v and the wandering exponent of the 
aperiodic sequence u [?,?]. 

In the case of a marginal perturbation, $ = 0, some exact results have been obtained for 
the Ising model. Critical exponents are found to vary continuously with the modulation 
amplitude [1,2,5-8] and, when the critical coupling is shifted, the bulk critical point 
becomes strongly anisotropic, i. e. the exponents of the correlation length, parallel 
and perpendicular to the layers, are different and their ratio varies continuously with 
the amplitude of the modulation [?]. 
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For some aperiodic sequences, the defect density vanishes in the thermodynamic limit, 
there is no shift in the critical coupling and the bulk critical properties remain unchanged. 
The Fredholm sequence [?], which belongs to this class and leads to a marginal perturbation 
for the layered Ising model, has been recently studied on a semi-infinite system [?]. 
Continuously varying surface exponents were obtained and it was shown that this type of 
aperiodic perturbation, which happens to be very regular, can be considered as a discrete 
realisation of the Hilhorst-van Leeuwen (HvL) model [?], i. e. a semi-infinite Ising model 
with smoothly varying couplings. 

In the present paper, we study the critical properties of a radial aperiodic perturbation 
in the two-dimensional Ising model. Instead of the parallel line defects of the layered 
system, we consider concentric circular defects distributed according to the Fredholm 
sequence. This type of perturbation is closely related to the radial HvL defect [11-13]. 

Our main motivation is to check the validity of the gap-exponent relation [?, ?] which 
follows from the transformation of the critical correlation functions under the conformal 
mapping of the original system onto a strip [?]. Since conformal transformations cannot 
be used on strongly anisotropic systems [?] , such a relation is excluded in the case of bulk 
aperiodic perturbations. But it is known to be satisfied in some marginally inhomogeneous 
systems at the critical point, provided the inhomogeneity is properly transformed, as shown 
for a narrow line defect [?, ?] and later for extended defects of the HvL type [?,?,?,?]. 

The structure of the paper is the following. In section 2, we present the model and 
recall the properties of the Fredholm sequence. In section 3, the local magnetization is 
studied by the corner transfer matrix method of Peschel and Wunderling [?] and the local 
critical exponent is compared to the value for the HvL defect. In section 4, the validity of 
the gap-exponent relation is discussed and our conclusions are given in section 5. 



2. The radial Fredholm perturbation 

We consider, in the (p, #)-plane, a two-dimensional Ising model with Hamiltonian 



in the continuum limit. Here H c is the critical Hamiltion of the unperturbed system, e 
is the energy density operator and the energy perturbation, with an amplitude g per unit 
length, is located on circles with radii p p = m p . 

The contribution of negative values of p to the perturbation is irrelevant since it 
renormalizes to a point defect with a finite amplitude. The local critical behaviour, which 
is governed by the long distance behaviour of the perturbation, remains unaffected if one 
considers circles with radii distributed according to the generalized Fredholm sequence [?] , 
i. e. with p p = m p + 1, p>0 and integer values of m>l . 

This aperiodic sequence, which is the characteristic sequence of the powers of m, follows 
from the substitution on the three letters A, B and C: 




= m 



v 



(2.1) 
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A -> S(A) = ABC C ■■■ C 
B^S(B) = B C C C ■■■ C 

C -> 5(C) = cccc ■■■ c 



B^S(B) = B C C C ■■■ C ^ ^ 



m 



With words of length to = 2, one recovers the usual Fredholm substitution [?]. Starting 
the substitution with A and associating to the k-th letter in the sequence the digit f k = 
for A and C, and /"^ = 1 for £?, after n iterations and with to = 2, one obtains: 



n = A 

n = 1 AS 

n = 2 A B B C 

n = 3 ABBCBCCC 

n = 4 ASSCSCCCSCCCCCCC 

/ fc 0110100010000000 



(2.3) 



which gives f k = 1 when = 2 P + 1. 

The substitution matrix has eigenvalues m, 1, 1 [?], so that the wandering exponent u> 
vanishes, leading to a marginal layered perturbation according to (??). It is easy to verify 
that the perturbation is also marginal with a radial defect. The total perturbation inside 
a circle with radius R ^> 1 is obtained by summing over the contributions of the circles up 
to Pmax ~ In .R/ In TO, 



'max 

Pmax 



E - S < 2 ' 4 ' 

and the average perturbation per unit surface at a length scale R is given by 

g(R) - | • (2-5) 

Under rescaling by a factor b = R/R', this thermal perturbation transforms according to 

L = »,.L (2 . 6 ) 

so that, with v = 1 in the two-dimensional Ising model, g is scale invariant, i. e. the 
radial aperiodic perturbation is marginal. Furthermore, according to (??), the average 
perturbation density vanishes in the thermodynamic limit, leaving the bulk critical point 
unchanged. 



3. Corner transfer matrix study of the local magnetization 



In this section, the critical behaviour of the local magnetization at the center of the 
defect will be deduced from a finite-size scaling analysis at the bulk critical point. This 



Figure 2. Through a rescaling of the lattice parameters, one restores isotropy near the 
critical point. The sector of figure f (dark triangle) has now a vanishing opening angle 
9. Covering the plane with such sectors, one generate the radial Fredholm defect (grey 
circles). 

can be achieved using the corner transfer matrix method of Peschel and Wunderling [?] 
which allows a study of rotation symmetric defects, while working on a lattice. 

We first consider the sector of a lattice Ising model shown in figure 1. With the same 
lattice parameter a± = a 2 = a in both directions, the opening angle is 9= ir/2. There are 
N + 1 horizontal rows with fixed boundary conditions on the last one. The interactions are 
between first-neighbour spins, constant and equal to K 1 in the horizontal direction, while 
they are given by K 2 (k) and depend on the row index, k = 0, N, in the vertical direction. 
This dependence follows the generalized Fredholm sequence, so that the vertical couplings 
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take the form 

K 2 (k) = K 2 r fk+1 , (3.1) 

where r is the modulation factor and the / fc s are the digits defined in equation (??). 

Let us take the extreme anisotropic limit, K 2 — > and K 1 — > oo, while keeping the 
ratio A = K 2 /K* fixed {K* = — \ lntanhi^ — > is a dual coupling). On the square lattice, 
the correlation lengths become different in the two directions with £ 2 /£i — <C 1 near 
the critical point. In order to recover an isotropic system which properly describes the 
continuum problem of the last section, one has to rescale the lattice parameters such that 
£i a i = ^2°2' which gives a 1 = 2K* — ^ if one takes a 2 for the unit lentgth. The opening 
angle of the sector is now reduced to 9 = 4K^ and the number of sectors needed to cover 
the plane, 

goes to infinity in the extreme anisotropic limit. 

In this way, as shown in figure 2, the rotation symmetry of the original problem is 
restored and the system becomes continuous along the defects. The perturbation per unit 
length is (r — l)K 2 /a 1 which allows us to identify the continuum parameter 

g=\{r-l) (3.3) 

at the critical point A c = 1. 

In the extreme anisotropic limit, the Baxter corner transfer matrix T [?] takes the 
simple form 

T = exp(-±eH) (3.4) 

where 9 is the infinitesimal opening angle of the sector and 7i is the Hamiltonian of the 
inhomogeneous quantum Ising chain: 

JV-l JV-l 

« = ~\ E 2k °l - \ E ( 2k + !) A ( fc ) °t°t+i ' A ( fc ) = Xrfk+1 ■ ( 3 - 5 ) 

k=l k=0 

The as are Pauli spin operators and the coupling X(k) has the same aperiodic modulation 
as K 2 in equation (??). 

The magnetization of the central spin, given by 

_ TrK^T") 

°~ Tr (T n ) ' [ * 

can be reexpressed in terms of the Fermi operators which diagonalize the corner transfer 
matrix [?, ?]. In this way, one obtains the local magnetization as a product, 

N 

m = Y\ tanh (^Tre^) , (3.7) 
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Figure 3. Variation of the local magnetization exponent (3i with the modulation factor r 
for the generalized Fredholm sequence with m = 2 and m = 3. The lines correspond to the 
conjectured expression of equation (3.9) and the points are finite-size scaling estimates. 



Table 1. Extrapolated finite-size estimates for the local magnetization exponent pi and 
the shift A contributed by the localized mode when r < 1, as functions of the modulation 
factor r for the Fredholm radial defect with m = 2 and m = 3. The figures in brackets 
give the estimated uncertainty in the last digit. Each second column gives the expected 
values of equation (3.9). 



r 








m = 2 








m = 3 






Pi 




A 




Pi 




A 




0.3 


1 


7373 


1.7370 


1.73697(6) 


1.73696559 


1.0973 


1.0963 


1.0959031(5) 


1.09590327 


0.4 


1 


3221 


1.3220 


1.32192811(4) 


1.32192809 


0.8373 


0.8356 


0.834043(1) 


0.83404377 


0.5 


1 


0010 


1.0006 


0.9999999(1) 


1 


0.6382 


0.6362 


0.630929(2) 


0.63092975 


0.6 





7405 


0.7398 


0.736964(3) 


0.73696559 


0.4801 


0.4780 


0.46493(6) 


0.46497352 


0.7 





5254 


0.5245 


0.51462(7) 


0.51457317 


0.3538 


0.3522 


0.326(2) 


0.32465952 


0.8 





3505 


0.3498 


0.32191(5) 


0.32192809 


0.2546 


0.2537 


0.211(6) 


0.20311401 


0.9 





2168 


0.2166 


0.1517(3) 


0.15200309 


0.1795 


0.1793 


0.101(4) 


0.09590327 


1.0 





1250 


0.125 






0.1250 


0.125 






1.1 





0692 


0.0690 






0.0870 


0.0868 






1.2 





0382 


0.0377 






0.0611 


0.0605 






1.3 





0215 


0.0208 






0.0436 


0.0425 







where the e^s are the N nonvanishing diagonal fermion excitations of the quantum chain|| 
(see reference [?] or appendix B in reference [?] for details) . 



|| The lowest excitation eo vanishes due to the fixed boundary conditions. 
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The excitations of the Hamiltonian (??) at the critical point A c = 1 have been obtained 
numerically on chains with N = m p + 1 spins with p = 4 to 16 for m = 2 and p = 3 to 9 
for m = 3. The modulation factor r has been varied from 0.3 to 1.3 with steps of 0.1. The 
critical value of m on a finite system vanishes as N~ Xl , where the scaling dimension of the 
local magnetization, x t , is continuously varying in the present case. This exponent is also 
equal to the local magnetization exponent f3 l since v = 1 in the 2d Ising model. The finite- 
size estimates for f3 l obtained from sequence extrapolations using the BST algorithm [?] 
are shown in figure 3. 

The layered Fredholm modulation on a semi-infinite system is known to lead to the 
same critical behaviour as the HvL model with couplings K 2 (k) — K 2 (l + a/k) if one 
makes the correspondence [?] 



a 



lnr 
lnm 



(3.8) 



Assuming the same relation for the radial defect and using the analytical result of 
reference [?] for the HvL model, we are led to the following conjecture for the Fredholm 
exponent^ 



A = 2 



lnr 
lnm 



du 



sinh 2 u 



sinh (2n I j 2 ^- 1 cosh it) 

V lnm I / 





lnr 


(* 


lnm 



(3.9) 



Here the quantity in brackets is a shift A, contributed by a localized mode, which has 
to be added when r < 1. Our numerical estimates for f3 l are in reasonable agreement 
with this expression as shown in table 1. The uncertainties on estimated by comparing 
different extrapolated values given by the BST algorithm, are much smaller than the 
actual deviations from the conjectured values. Such a behaviour is known to occur when 
logarithmic corrections to scaling are present [?]. The correspondence (??) is strongly 
supported by the extrapolated values of the shift A for which the agreement with the 
expected values is excellent. 



4. Conformal mapping and gap-exponent relation 

Let us now consider the transformation of the marginal perturbation in the continuum 
limit, equation (??), under the logarithmic conformal mapping [?] 

w = — In z , w = u + iv , z = pe li} , (4.1) 

which transforms the whole plane into a strip (— oo < u < 4oo, < v < L) with periodic 
boundary conditions in the r-direction. The local dilatation factor is b(z) — \ dwfdz\~ 1 = 



% The amplitude a is half the one used in reference [?]. 



Figure 4. The logarithmic conformal mapping transforms the radial aperiodic Fredholm 
defect in the plane into a periodic defect on the strip. The period L m is proportional to 
the strip width L. 



2irp/L and the amplitude of the thermal perturbation, t(z) = g^2 p 8(p — p p ), is changed 
into: 



' n— — no ^ 



2nu 



- mP 



2n 



1/1/ +00 



exp 



2tvu 



(4.2) 



1 

1 - - 

v 



exp 



p mm 



27TW 



With the Ising model, v = 1, the perturbation is marginal, and the first part of the -in- 
dependence is eliminated so that: 



t(w) = g 



2n 

T 

+00 



+00 



S<l-exp 

=—00 ^ 

E. / Llnm\ 
6 { U - p — 



2tt 

T 



u — p ■ 



Llnm 



2tv 



(4.3) 



As shown in figure 4, the perturbation now consists of straight line defects, the distance 
between two successive lines, 



L m = Llnm/2n , 



(4.4) 



being proportional to the width of the strip L. The discrete dilatation invariance of the 
radial defect on the plane has been transformed into a discrete translation invariance along 
the strip. 
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In the case of the HvL radial defect, the continuous dilatation invariance of the 
perturbation, t(z) = g/p on the plane, leads to a constant deviation from the critical 
coupling t(w) =g(27c/L) on the strip [?, ?, ?]. If £^ denotes the correlation length associated 
with the spin-spin or the energy-energy correlations on the off-critical isotropic strip, it 
satisfies the following finite-size scaling relation [?] 

^ 1 (t,L) = L- 1 X (j) (cL 1 /»t), (4.5) 

where the gap scaling function X^{t) is universal [?], t is the deviation from the bulk 
critical temperature and c is a nonuniversal constant. The defect scaling dimension follows 
from the gap-exponent relation [14-16] with: 

xt = ^ 1 (t,L) = ^-X^(cL^t). (4.6) 

It is continuously varying with the defect amplitude g, as expected for a marginal 
perturbation, but it varies in a universal way: different models belonging to the same 
class of universality will show the same dependence. 

With the Fredholm radial defect, which gives a periodic layered perturbation in the 
strip geometry, one could directly deduce the gap (inverse correlation length) from the 
appropriate transfer matrix on the strip and use (??). Alternatively, one may estimate the 
deviation t from the bulk critical temperature in (??) which is associated with the periodic 
defect. A measure of this deviation is provided by the shift in the critical temperature 
induced by the same perturbation (i. e. parrallel line defects with a fixed distance L m ) on 
the infinite plane. 

In order to avoid the calculation of the nonuniversal constant c, we shall proceed via 
a comparison to the HvL problem. We use the same extreme anisotropic limit as for the 
corner transfer matrix, with vertical couplings K 1 ^oo, horizontal couplings K 2 — >0 and 
vertical ladder defects corresponding to modified couplings 

K' 2 = rK 2 (4.7) 

as shown in figure 4. As above, a length rescaling (a 2 = 1, a x = 2K^) is understood 
in order to restore isotropy. For the radial HvL problem the perturbed interactions 
K 2 (k) =K 2 (1 + a/k) lead to an homogeneous shift of the horizontal couplings with 

K' 2 = K 2 (l + aj?) . (4.8) 

The critical temperature of a periodic layered anisotropic system follows from the 
relation [?] 

2> 2 *0-) = £*itf)' (4.9) 

3 3 

where K 2 -bonds are perpendicular to the layers and the sums are over a period. In the 
case of the Fredholm problem, this leads to 

(L m - 1)K* + K' 2 * = L m K x , (4.10) 
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where, in the extreme anisotropic limit, equation (??) leads to 



K' 2 * = K* - \\nr . 
Putting (??) into (??) gives the criticality condition: 
lnr 



Kt 



2L 



= K l . 



(4.11) 



(4.12) 



The corresponding relation for the HvL problem is obtained using (??) with L m = l. Since 
the modified coupling in (??) corresponds to r = l + a(2n/L) : with L» 1 equation (??) 
gives: 



K'* « K* 2 - a- 
and finally the criticality condition reads: 



K* 



a I = K 1 . 



(4.13) 



(4.14) 



Comparing equations (??) and (??) and using (??) leads to the conjectured correspondence 
of equation (??). 

Besides the local magnetization exponent of equation (??), the gap-exponent relation 
also gives the local energy exponent of the radial Fredholm defect as [?] 



x. 



1 + 2 



lnr 
mm 



+ o 



lnr 
lnm 



(4.15) 



where, due to the self-duality of the Ising model, only even powers of lnr/ lnm enter the 



expansion 



5. Conclusion 

The close relationship between the discrete aperiodic Fredholm defect and the 
continuous HvL inhomogeneity in the two-dimensional Ising model has been established 
in a study of the surface critical behaviour of the layered system in [?]. It has been also 
recently verified for the bulk layered Fredholm defect [?]. 

In the present paper the validity of this connection -which is summarized in 
equation (??)- is further extended to radially symmetric defects. The inhomogeneity in this 
case corresponds to an infinite sequence of concentric circles with exponentially increasing 
radii. The separation of distant circles becoming very large, the density of defect bonds 
goes to zero, thus the bulk critical behaviour of the system remains unchanged. However 
an infinite sequence of defect circles modifies the local critical behaviour at the centre of 
the defect. 
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Our second observation, concerning the validity of the gap-exponent relation for the 
Fredholm defect, is somewhat unexpected. It was known till now [?] that some aspects 
of conformal invariance, including the gap-exponent relation, are still satisfied for some 
marginally inhomogeneous systems, which either contain a finite number of defect lines 
or display a smooth variation. In the Fredholm problem the number of defect circles is 
infinite, furthermore the perturbation changes the continuous dilatation symmetry into a 
discrete one. 

A similar analysis of the layered Fredholm defect problem [?, ?] reveals that the gap- 
exponent relation stays valid in this case, too. Now the transformed perturbation on the 
strip is still periodic with period L m as given in equation (??), however its shape is much 
more complicated than for the radial defect. Taking the anisotropic limit and considering 
the product of L m successive transfer matrices one can establish the same relation with 
the HvL model as in (??). These observations lead us to conjecture the validity of the gap- 
exponent relation for any marginal perturbation which does not modify the bulk critical 
behaviour of the system. 

If the marginal perturbation extends over the volume of the system with a nonvanishing 
density, the gap-exponent relation is generally no longer valid. As mentioned in the 
introduction, this type of marginal aperiodic perturbations lead to strongly anisotropic 
systems [?] which cannot be transformed using conformal techniques. 
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